Evolutionary Applications

Restoration and management of natural populations often assume that local genotypes are best suited for transplantation to their local environment. Prioritizing a single local genotype, however, contrasts with the framework of maximizing intraspecific diversity to increase population resilience to environmental change. Local populations may also become maladapted to a rapidly changing environment, motivating alternative frameworks that instead minimize environmental distance between source and transplantation sites. Here, we tested the predictive power of the local is best, maximize intraspecific diversity, and minimize environmental distance frameworks on the survival and growth of Eastern oyster (Crassostrea virginica) genotypes in field common gardens that differed in salinity and disease pressure. Although a genome scan revealed patterns of adaptation to disease, heat stress, and salinity among source populations, we did not find strong support for the local is best framework: geographically distant southern genotypes performed comparably to local selection lines and a local wild population. Higher genetic diversity within monocultures was associated with higher survival, yet highly diverse polycultures survived at lower rates than the best-performing monocultures, providing mixed support for the maximize intraspecific diversity framework. Temperature and salinity of the environments-of-origin of parents predicted the survival of their offspring in common gardens, but the relationship between survival and environmental distance was context-dependent, leading to mixed support for the minimize environmental distance framework. Together, these results demonstrate that no single framework reliably predicted transplantation success, suggesting that effective management strategies may need to integrate genomic and environmental lines of evidence to guide genotype selection.

Trait-selective harvesting by fisheries can impose strong selective pressures on fish populations, driving changes in life history traits affecting fisheries productivity and ecosystem functioning. While the genetic consequences of harvesting have been extensively studied, the extent to which phenotypic variation reflects genomic evolution versus environmentally-induced plasticity remains unclear. Epigenetic mechanisms, such as DNA methylation, may mediate between these processes, serving as a rapid and reversible response to the selective pressures imposed by harvesting. In this study, we implemented an improved laboratory and bioinformatics protocol, epiGBS3, to examine genomic variation and DNA methylation patterns in the marine fish Xyrichtys novacula. The study spanned three replicated geographical areas each comprising two adjacent locations: an intensively exploited fishery and a no-take Marine Protected Area (ntMPA). A nested analysis design across the three areas revealed strong gene flow and no evidence of genetic structure. Nevertheless, nucleotide diversity was significantly reduced in fisheries relative to ntMPAs. We also found that DNA methylation levels differed between protected and exploited sites after controlling for age, suggesting that fishing may influence epigenetic changes independently of fisheries-induced age-truncation effects. This represents one of the first lines of evidence that fisheries can potentially shape epigenetic variation, supporting DNA methylation as contributor to local adaptation under high gene flow and strong anthropogenic selection.

Selective breeding in aquaculture is necessary to establish food security and meet demand for sustainably produced protein. An informed selective breeding program requires understanding how population structure, environmental adaptation, and human activities shape natural genetic variation in wild conspecifics. Unfortunately, wild variation remains poorly characterized for many commercially important aquaculture species. Here, we conduct the first range-wide study of genomic population structure for the eastern oyster (Crassostrea virginica) across thousands of miles (Texas, USA to Eastern Canada) using a 200K SNP array. We integrate population structure analyses, genotype-environmental associations, and structural variant detection to identify adaptive loci and quantify human-mediated genetic impacts. Our data confirms two ancestral clusters with a phylogeographic break between the Gulf and Atlantic (FST = 0.06) and highlights patterns of substructure within each region. We find evidence of unexpected patterns of genomic variation in two locations: evidence of Gulf ancestry in a mid-Atlantic estuary (Chesapeake Bay), and evidence of Atlantic ancestry in a Gulf estuary (Apalachicola Bay). While we cannot definitively determine the causes of these unexpected patterns, we show that they are consistent with direct and indirect human impacts in these estuaries. Genotype-environment association analyses with in situ temperature and salinity measurements were used to identify putatively adaptive loci, including SNPs within large structural variants (>1Mb). Our results identified genomic targets for aquaculture breeding programs aimed at climate resilience, reveal complex patterns of human impacts in managed systems, and demonstrate how seascape genomics can be used to improve aquaculture outcomes.

Successful establishment of a species in a new range is a useful way to understand the impact of demography and selection on the evolution of globally distributed species. In particular, introductions influence genetic diversity and population structure in the introduced range in unpredictable ways. Additionally, introgressive hybridization is often associated with successful establishment in new ranges. In this study, we explore the impact of introgressive hybridization on the polyploid Capsella bursa-pastoris in the New York City metropolitan area. We find Capsella bursa-pastoris in the New York City metropolitan area likely originated from multiple introductions from northern Eurasia, and that populations across the New York City metropolitan area are generally panmictic. As with Capsella bursa-pastoris in Eurasia, we discover evidence of introgression from the diploid Capsella rubella in this population. By evaluating ancestry in regions across the genome, we find introgressed regions are rich in gene content and contribute to genetic diversity in this population. These results suggest that introgressive hybridization before introductions may buffer species from the negative effects of population bottlenecks and allow for successful establishment.

Foodborne non-typhoidal Salmonella remains a major public health concern, yet routine surveillance recovers large numbers of isolates from food that are not associated with human illness. Studies have shown foodborne isolates can be genetically linked to clinical cases, highlighting a critical challenge for risk assessment and outbreak prioritisation. This study aimed to determine whether genomic markers can distinguish foodborne Salmonella strains with an increased likelihood of causing infection. Whole-genome sequencing data from over 900 Salmonella isolates recovered from food and the environment through UK Health Security Agency surveillance were analysed using hierarchical clustering to define genetically related groups. These clusters were expanded using the global EnteroBase database to provide broader epidemiological context. Genome-wide association analyses identified genetic markers associated with clusters containing clinical isolates, including phage-associated regions. A highly conserved 7 kb marker identified in S. Agona demonstrated strong predictive performance at a global scale, with high sensitivity and specificity for infection-associated lineages and strict serovar restriction. Comparative genomic analysis revealed that all markers localised to a shared chromosomal hotspot corresponding to a prophage integration site. The 7 kb risk-associated marker formed part of a larger prophage closely related to the well-characterised S. Typhimurium Fels-2 phage, which encodes a DNA invertase linked to phase variation, a mechanism known to promote phenotypic heterogeneity and host adaptation. As these S. Agona isolates are monophasic, our findings indicate that our genome-wide association approach has rediscovered this DNA invertase known to contribute to infection risk but in a different serovar via an alternative regulatory mechanism. Overall, this work demonstrates the potential to move beyond treating all foodborne Salmonella isolates as equivalent hazards, towards a genomics-informed framework for risk stratification. This approach provides a foundation for improved risk-based decision-making, enhance outbreak investigations and enable earlier prioritisation of public health responses during Salmonella surveillance and control. Author summaryFoodborne Salmonella infections remain a major public health concern, but not all strains pose the same risk to human health. Here we investigated whether genetic differences could explain why some foodborne strains are more likely to cause human infection. We analysed over 900 genomes from food and environmental sources, grouping closely related strains before placing them in a global context using EnteroBase. By combining pangenome and genome-wide association analyses, we identified distinct lineages within several serovars that differed in their association with human cases. In Salmonella Agona, all clinical isolates belonged to a single lineage carrying a highly conserved 7 kb marker that was absent from low-risk strains. This marker demonstrated strong sensitivity and specificity across global datasets and was located within a prophage closely related to the well-characterised Fels-2 phage. This region encodes a DNA invertase previously linked to phase variation, a mechanism that promotes bacterial adaptability. Our findings indicate that infection risk can be structured at the lineage level and influenced by mobile genomic elements, particularly prophages, that enhance environmental persistence and host adaptation. This work advances genomic surveillance from retrospective linkage towards mechanistic and predictive risk assessment, with direct relevance for supporting risk-based decision-making during outbreak investigations.

Sea star wasting disease has caused widespread mortality in the kelp forest predator, the sunflower sea star (Pycnopodia helianthoides). Wild populations have declined by up to 99% in parts of their native range along the western North American coast. In response, a multi-institutional conservation breeding and rearing program has been initiated to support future reintroduction efforts for the species. We split a full-sibling cohort across four larval density treatments (1 larva/ml, 2 larvae/ml, 5 larvae/ml, and 15-20 larvae/ml) to assess the effects on larval settlement, juvenile survival, and juvenile fitness at 12 months old. Stars raised in the highest density treatment displayed a lower settlement rate and were significantly smaller than the other density groups at 12 months old, but showed no significant difference in flip time, a measure of fitness. Additionally, measurements of diameter, weight, and arm count across modern and historical juvenile and adult stars indicate that P. helianthoides experience exponential weight gain as they grow in length, with corresponding asymptotic growth in arm count. These findings will inform best practices for the aquarium propagation of P. helianthoides and will contribute to broader efforts aimed at reestablishing populations in the wild.

The plant-beneficial bacterium Pseudomonas protegens CHA0 (CHA0) is widely studied for the biological control of soil-borne plant diseases. Beyond its root-colonising capabilities, CHA0 can also infect and kill insect larvae and thus exhibits a multi-host lifestyle shared with other plant- and insect-colonising bacteria. To better understand the robustness of this multi-host lifestyle, we subjected CHA0 to ten consecutive passages through larvae of the pest insect Plutella xylostella via repeated cycles of insect colonisation and killing forcing it into an insect-only lifestyle. Overall, serial passaging did not result in consistent changes in insect killing speed, larval or root colonisation, plant protection efficiency, microbial antagonism or in vitro growth. This suggests that its multi-host lifestyle was conserved following serial passage. Nonetheless, a few independently passaged lines showed an increase in larval killing speed, which in one case might be linked to choline uptake. To disentangle changes specific to the insect host from those arising due to the experimental system itself, we conducted parallel serial passages through the same system while omitting the insect host. In some of these lines, exposure to the background of the system led to changes in microbial antagonism and in in vitro growth, which likely are associated with mutations in regions encoding for regulatory systems. Our findings indicate that P. protegens CHA0 remains phenotypically stable in complex environments such as an insect host, suggesting that the multi-host lifestyle might also be conserved when applied in the field and supporting CHA0s potential for reliable biocontrol performance against both plant diseases and insect pests. Author summaryControlling insect pests with living organisms, known as biological control, offers an environmentally friendly alternative to chemical pesticides. The plant-beneficial bacterium Pseudomonas protegens CHA0 is a promising biocontrol candidate that not only colonizes plant roots but also infects and kills certain insect larvae. This ability to colonize different hosts appears to be a conserved trait also observed in other bacteria. To better understand the robustness of this multi-host lifestyle, we repeatedly exposed CHA0 to larvae of the insect pest Plutella xylostella and assessed the resulting physiological and genetic changes. Surprisingly, after ten cycles, CHA0 largely retained its insect-killing and plant-protective traits. Although a few populations showed minor changes, including slightly faster insect killing and traits associated with aspects of the experimental system, these changes were limited in scope. Overall, our findings suggest that P. protegens CHA0 does not change rapidly in complex environments such as an insect host, supporting its potential for reliable biocontrol performance in the field.

Introgressive hybridization between wild and domestic animals is a widespread phenomenon with important implications for genetic diversity, local adaptation, and conservation management. The causes and consequences of this process are poorly understood. In Australia, hybridization between dingoes and domestic dogs presents a dual conservation challenge, threatening the genetic integrity of dingoes while allowing potential adaptive introgression. To investigate the environmental drivers of this process, we analyzed high-density SNP array data in 390 dingoes and 396 domestic dogs. Dingo populations showed regional genetic structure and were clearly differentiated from domestic dogs. Using local ancestry inference and genome-environment association analyses, we found low levels of dog introgression in dingoes from remote areas in Central and Western Australia, and moderate levels in Eastern and Southern populations. Climatic variables (maximum temperature of the warmest month, mean temperature of the driest quarter) and the Human Footprint Index (reflecting density of human populations and environmental modifications) were significant predictors of introgression. We identified four genomic regions with overrepresented dog ancestry, including a large introgressed block on chromosome 27, which contained an olfactory receptor gene showing signatures of positive selection, suggesting adaptive introgression. In addition, a chromosomal inversion previously described in dogs and absent in dingoes was initially identified as an introgressed block. We also detected eight genomic regions nearly free of dog ancestry, suggesting purifying selection against maladaptive variants. Together, these results highlight the complex interplay between introgression, human influence, and local adaptation in dingoes, offering valuable insights for conserving the evolutionary potential of this apex predator in increasingly modified landscapes.

Kelp forests are widely distributed along temperate and polar coastlines worldwide and are among the worlds most productive and diverse marine ecosystems. Yet, due in part to ocean warming, they are declining and even disappearing in many parts of the world. While genomic tools can identify local adaptation and predict species responses to global change, these predictions have rarely been validated in the field, hampering their widespread use in conservation practice. Here, we applied a seascape genomics approach to investigate environmental adaptation in the two main canopy-forming species of the Northeast Pacific, Macrocystis tenuifolia and Nereocystis luetkeana. We leveraged whole-genome sequences of 598 individuals across 94 sites along the British Columbia and Washington coasts, together with 37 environmental variables. Both species showed genomic signatures of local adaptation, with distinct environmental drivers shaping adaptation in each species despite their co-occurrence across much of the studied area. Using gradient forests, we modelled the genetic turnover across environmental gradients and predicted populations vulnerability (genomic offset) under projected environmental conditions. Genomic offsets differed greatly among regions and were positively correlated with kelp declines observed to date, especially in Macrocystis, validating the link between genomic models and outcomes in the field and allowing us to translate genomic predictions into an ecologically meaningful metric: the risk of extirpation under global change. Our models predict that assisted migration could significantly attenuate kelps vulnerability to global change. Across environmentally heterogenous coastlines, short-distance migration can often substantially reduce future genomic-environmental mismatches, but in many cases, long-distance migration would be most beneficial. Our results highlight the potential of seascape genomics to predict vulnerability of populations to global change. Importantly, the validated link between our genomic models and ecological outcomes allows quantification of climate-driven extirpation risk and can inform conservation strategies to improve the resilience and sustainable management of these vulnerable ecosystems.

BackgroundThe virulence-persistence trade-off is considered a fundamental organizing principle of Listeria monocytogenes population biology: hypervirulent clonal complexes dominate clinical cases but are rarely found in processing environments, while hypovirulent lineages dominate industrial niches but are underrepresented in severe disease. However, whether this dichotomy operates as an absolute paradigm has not been quantitatively evaluated at the population scale. Here we develop a multi-dimensional genomic scoring approach that simultaneously quantifies virulence potential (V), environmental persistence capacity (P), clonal epidemiological context (C), and antimicrobial resistance (R) across 903 genomes from four independent datasets spanning five countries, and apply it to test the universality of the trade-off and to characterize the ecological strategies of L. monocytogenes at the population level. MethodsThe scoring approach integrates four components into a composite 0-100 score through empirically calibrated weights (V: 30%, P: 40%, C: 20%, R: 10%). Validation employed 903 L. monocytogenes genomes from four public BioProjects: longitudinal industrial surveillance in Norway (Fagerlund et al. 2022, n = 513, PRJNA689484), retail environments in the United States (Stasiewicz et al. 2015, n = 191, PRJNA245909), clinical-environmental context in China (Wang et al. 2021, n = 151, PRJNA759341), and meat processing in Poland (Kurpas et al. 2020, n = 48, PRJNA629756). ResultsThe composite score achieved excellent discriminatory performance for identifying persistent clones (AUC = 0.933; 95% CI: 0.910-0.954) with perfect specificity (1.000; zero false positives). The inverse V-P correlation was statistically significant across all four datasets (Spearman {rho} from -0.144 to -0.713; p < 0.01), providing the first cross-dataset quantitative confirmation of the trade-off. However, simultaneous evaluation of V-P profiles at the population scale revealed that the species does not conform to a binary dichotomy but rather exhibits three quantitatively distinguishable ecological strategies, for which we propose a functional trophic taxonomy: nosotrophic lineages (22.7%; V > 65, P < 35), specialized in the pathogenic niche; saprotrophic lineages (5.8%; V < 30, P > 45), with irreversible virulence attenuation and industrial specialization; and, as the central finding, amphitrophic lineages (39.1%; V [&ge;] 35, P [&ge;] 40), which simultaneously retain functional inlA and stress tolerance determinants (SSI-1) without detectable genomic sacrifice. The three strategies differed significantly (Kruskal-Wallis H = 138.7; p = 7.6 x 10-3{superscript 1}). The correspondence between trophic strategy and CC was predominant but not absolute, demonstrating that this phenotypic classification captures intra-CC functional heterogeneity inaccessible through conventional typing. Furthermore, comparison between genome-based and surveillance-informed classifications revealed that 60 hypervirulent isolates (CC1/CC14), genetically classified as nosotrophic, persisted for up to 8 years in industrial facilities despite lacking any recognized persistence markers -- indicating that their prolonged survival reflects environmental opportunity rather than intrinsic genomic adaptation. ConclusionsMulti-dimensional genomic profiling reveals that the virulence-persistence trade-off, while statistically robust, does not operate as an absolute paradigm. The amphitrophic strategy -- documented here for the first time as a quantitatively distinguishable category encompassing 39.1% of the analyzed population -- challenges the prevailing dichotomous model and identifies a previously unrecognized combined ecological niche. The ability to discriminate between genome-encoded persistence capacity and environmentally facilitated persistence provides a biological framework for understanding the ecological determinants of L. monocytogenes population dynamics in anthropogenic environments.

When invasive populations establish in regions far from their origin, they risk accumulating harmful mutations (genetic load) that limit population viability and subsequent spread. This may be exacerbated by the multiple, sequential bottlenecks experienced when invasions stem from a bridgehead population. However, populations may be able to purge genetic load when they can outbreed with other lineages from subsequent invasions. Here, we analyse global invasions of a species complex of persistently inbreeding ambrosia beetles, using genomic data (N=247) from invasive populations in Africa, North America and Australia, and from native populations in Asia. We focus particularly on one species of this complex (Euwallacea fornicatus) which poses a catastrophic threat to tree species worldwide and is rapidly expanding its global range. We uncover a single lineage of this species across South Africa, California and Western Australia, derived from an invasive bridgehead and containing almost no nuclear genetic variation. In South Africa we identify a second lineage that has repeatedly hybridised with the first lineage. Genetic patterns in the native range indicate that such opportunistic outbreeding may be common. Although purifying selection was evident in all lineages, native populations had fewer missense mutations than invasive populations, suggesting that opportunistic outbreeding may help purge fixed deleterious mutations when local lineage diversity is high. These findings show how inbreeding depression can affect populations even where inbreeding is common, and they highlight the biosecurity threat posed by subsequent gene flow into invasive populations.

The western honey bee (Apis mellifera) is an essential pollinator facing unprecedented threats from pesticide exposure. While pesticide resistance evolution is well documented in agricultural pests, our understanding of genetic variation in honey bee detoxification systems remains limited. This represents a missed opportunity, as harnessing naturally occurring detoxification diversity could provide new avenues for pollinator protection. Cytochrome P450 monooxygenases (CYPs), which are central to xenobiotic metabolism, offer a promising starting point. Here, we present the first comprehensive analysis of CYP genetic diversity in A. mellifera. We analysed the CYPome of 1,467 individuals representing 18 A. mellifera subspecies from 25 countries and identified 5,756 single-nucleotide polymorphisms (SNPs) in 46 CYP genes. Imputed McDonald-Kreitman testing revealed that 56% of non-synonymous CYP substitutions were driven by positive selection. Of the 1,302 haplotypes identified, 84% resided in CYP3, concentrated in the CYP9 and CYP6AS subfamilies implicated in xenobiotic detoxification. Population-level analysis of nucleotide diversity, Tajimas D selection signatures, FST-based differentiation, and McDonald-Kreitman testing pointed to CYP3 clan genes as the primary locus of adaptive variation. This work provides the first step toward building a comprehensive pharmacogenomic resource for honey bees, enabling the prediction of population-specific pesticide vulnerabilities and leveraging naturally occurring detoxification variants to enhance pollinator resilience - a critical step toward sustainable pollinator management.

Identifying spatial and temporal patterns of connectivity among populations is fundamental to marine ecology, evolutionary biology, and fisheries management. Yet, due to large population sizes and low genetic differentiation among populations, empirical quantification of population connectivity across a species entire range has not been achieved for an open-coast marine organism. Here, we leverage experimental transcriptomics to develop a genotyping-in-thousands by sequencing (GT-seq) panel to support assignment of recruits of the kelp forest gastropod, Kellets whelk (Kelletia kelletii), collected across the species biogeographic range. Over a three-year period, we identified high self-recruitment in the historical range (100%) and low self-recruitment in the expanded range (10.53 - 13.73%). Additionally, self-recruitment within the expanded range generally increased with recruit age, from 27.14% at 0.93 years to 43.40% at 1.93 years, indicating that the locally spawned individuals were more likely to survive to older ages than migrants from the historical range. Together, these results reveal limited self-recruitment in the expanded range and suggest that a post-settlement selective filter contributes to differential survival in a high gene flow marine system.

Nakaseomyces glabratus is a globally distributed opportunistic fungal pathogen. An ongoing discussion in studies of N. glabratus population structure has been whether genetic clusters are best defined using multilocus sequence typing (MLST) or short-read whole-genome sequencing (WGS). To assess the concordance between MLST- and WGS-based phylogenies, we analyzed a dataset of 548 N. glabratus WGS sequences from 12 countries. Clusters identified from WGS largely recapitulated the MLST-defined sequence type (ST) groups: fourteen WGS clusters were composed of a single MLST ST, and the remaining contained STs with very closely related MLST profiles. We thus propose a pragmatic naming convention, consistent with the system used in other microbial species, which specifies WGS cluster labels based on the primary ST. From the large WGS isolate dataset, we determined the prevalence of admixture and genomic variants. Interestingly, seven of the nine singleton isolates were admixed, in addition to 58 isolates from six different clusters. Aneuploidy was detected in 4% of isolates, most commonly in chrE, which contains ERG11, the gene encoding the enzyme targeted by azole antifungals. Aneuploid chromosomes did not exhibit elevated heterozygosity relative to the sequencing error rate, consistent with instability of extra chromosome copies. Copy number variants were found in 3% of the isolates; some of the CNVs co-occurred with aneuploidies, and were primarily identified on chrD, chrE, chrI, and chrM. Our findings demonstrate that deep splits between clusters preserve the utility of MLST ST designations for clade-level designation, yet underscore the utility of WGS for high-resolution genomic analyses. Article SummaryThere is an ongoing debate in studies on Nakaseomyces glabratus about whether traditional MLST analysis is sufficient to determine population structure, or whether the precision of whole genome sequencing (WGS) is necessary. We analyzed WGS data from 548 isolates from around the world. We found a very strong agreement between the two methods. We propose a hybrid naming system, where cluster names are based on the dominant MLST group. We used the WGS data to show that admixed isolates, and those with extra chromosomes or CNVs are rare (<7% of isolates in each class) and are distributed throughout the phylogeny.

Global networks of crop breeding programs leverage diverse germplasm, but diversity increases the complexity of maintaining stability in their elite genepools. To characterize genetic heterogeneity in breeding metapopulations and develop insights on how to manage it, we simulated the evolution of breeding populations on fitness landscapes. We revealed the geometric decrease in the average effect size of alleles segregating as standing variation that become fixed along an adaptive walk. We also demonstrated how independent adaptive walks of subpopulations are influenced by genetic drift, leading to cryptic genetic heterogeneity among elite genepools. This variation is released when elite lines derived from independent subpopulations are crossed, leading to segregation for 2-4X more major QTL in admixed families as in unadmixed families, and 2-4X more epistatic interactions. The emergent property of fitness epistasis for traits under stabilizing selection is well-understood in evolutionary genetics, but under-appreciated in crop quantitative genetics. To highlight the importance of this phenomenon, we constructed an empirical genotype-to-fitness landscape from the sorghum NAM, a global admixed prebreeding resource, demonstrating the utility of fitness landscapes for inferring genetic compatibilities within metapopulations. Our findings suggest that in breeding networks, strategies for effective germplasm exchange must account for epistasis in the oligogenic component of the genetic architecture of locally-adapted traits. Article summaryModern public sector crop improvement happens in networks of breeding programs that routinely exchange genetic information. Traditional models for understanding quantitative traits have limited predictiveness in situations with such genetic heterogeneity. This study uses breeding simulations and empirical data to show the utility of the fitness landscape framework for characterizing the genetic architecture of complex traits in breeding metapopulations. By simulating the evolution of breeding programs and integration into networks, it demonstrates how epistatic interactions between large-effect alleles are a fundamental property that must be accounted for when exchanging germplasm. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=102 SRC="FIGDIR/small/712732v1_ufig1.gif" ALT="Figure 1"> View larger version (25K): org.highwire.dtl.DTLVardef@1541326org.highwire.dtl.DTLVardef@b553a8org.highwire.dtl.DTLVardef@8758b4org.highwire.dtl.DTLVardef@1d0bdcd_HPS_FORMAT_FIGEXP M_FIG C_FIG

The ecological and evolutionary processes determining species range limits remain poorly understood. Ultimately, range limits depend on the species abilities to persist under heterogeneous conditions, by adaptive differentiation and phenotypic plasticity, including transgenerational effects. To investigate ecological differentiation and transgenerational effects in the clonal invasive knotweed, Reynoutria japonica, in Europe, we conducted a two-phase transplant experiment: plants sampled along the entire latitudinal gradient were planted in three sites located at the northern range margin, mid-range and near the southern range margin, and then re-transplanted among all three sites after two years. Biomass production and allocation were generally not associated with latitude of origin and previous growth at the same site did not promote performance. We therefore find no evidence that adaptive differentiation or transgenerational effects contribute to the wide distribution of R. japonica in Europe. However, at the northern site, with a 25% shorter season, knotweed plants invested much less biomass below-ground, and the pattern was further strengthened in plants that had grown in the northern site in the previous generation. Overwintering below-ground rhizomes are essential for survival and spread. We further explored limiting climate conditions in a species distribution model for the European range and found that mean annual temperature and temperature annual range are the main predictors of the European distribution of R. japonica. Taken together, our study suggests that low temperatures and associated short seasons may pose a limit to the broad environmental tolerance of R. japonica and restrict its northward spread by reducing below-ground biomass accumulation.

Malaria remains a major global health burden, with traditional control methods facing challenges such as insecticide resistance and high operational costs. Genetic biocontrol offers a promising alternative for mosquito population suppression, but its field efficacy would require assessment. This study evaluates the role that population genomic statistics can play in detecting decreases in population size in the context of a cluster randomized control trial (cRCT), investigating the response of nucleotide diversity ({pi}), Tajimas D, segregating sites, and linkage disequilibrium (LD) under both constant and seasonal demographic scenarios. We simulated 90% and 99% population declines with various degrees of between-cluster heterogeneity, and assessed the detection power of each statistic over time and number of clusters per arm. Results show that Tajimas D is highly sensitive and robust across crash severity, seasonality and heterogeneity scenarios. Segregating sites has similar power to Tajimas D when baseline data are available. We further estimated that cRCTs require approximately 3 to 5 villages per treatment arm to achieve adequate statistical power. These findings provide recommendations for genetic monitoring of vector control interventions in wild populations.

Population genetic theory predicts that natural selection will be more efficient in large than small populations because genetic drift reduces the efficiency of selection in small populations. Small populations adapting to new environments can also be expected to evolve higher recombination rates to facilitate adaptation as well as to dissociate and purge harmful mutations. We tested these hypotheses (1) by investigating differences in the strength of association between nucleotide diversity ({pi}) and recombination rate across the genomes of nine-spined sticklebacks (Pungitius pungitius) from four small freshwater (mean Ne {approx} 2 578) and four large marine (mean Ne = 86 742) populations, as well as (2) by comparing recombination rates between small and large populations using population specific linkage maps. We found the predicted positive correlation of{pi} with recombination rate from all but the smallest freshwater populations, suggesting prevalent linked selection even after accounting for variation in GC/CpG content, and gene density. Mean recombination rates did not differ between freshwater and marine populations, except that the smallest Ne freshwater population exhibited significantly elevated recombination rate. GWAS analyses suggested a polygenic basis for recombination rates. These results suggest an important role for linked selection in reducing{pi} in low recombination regions especially in large populations. Moreover, as predicted by theory, at least one of the small freshwater populations appears to have evolved a higher recombination rate than its marine ancestors.

Coral reefs face widespread declines under global change, yet the status of their genetic diversity remains largely unknown. Here, we compiled genome-wide DNA sequencing data for 2,520 individuals from 18 reef taxa--including corals, fish, sharks, oysters, shrimp, sea anemones, and manta rays. These data were used to assess spatiotemporal patterns of genetic distances across 173 reefs worldwide between 1998-2018. While we did not observe an overall temporal decline in genetic distances, within-reef distances showed negative temporal trends, potentially reflecting population-level diversity loss. These effects varied across species but did not show clear distinctions between taxa. We then used satellite-derived seascape variables to predict local effects on genetic distances across reefs globally. Negative effects were predicted for the Red Sea, Northern Caribbean, and Coral Triangle, while positive effects were found across the South Pacific. Key predictors included declining oxygen levels, increasing nitrate concentrations, and rising water temperatures--variables that can be tracked in real time via Earth observation (EO), enabling early warning for coral reef genetic diversity loss.

Bacteriophages are not only the most ubiquitous biological entity on earth, they also display remarkable genetic diversity across and within populations. While macrodiversity has been extensively studied, the drivers of microdiversity (intraspecies genetic diversity) remain poorly understood, particularly in relation to phage lifestyle. The distinguishing ability of temperate phages to integrate themselves into the host genome has an unknown influence on the microdiversity present. This difference in microdiversity could impact the adaptability of phages to (a)biotic factors. To identify a possible association between microdiversity and lifestyle, we analysed 12 existing viromics datasets focusing on soil bacteriophages, including 41 412 viral genomes in total. We found that phages predicted to be temperate consistently exhibit significantly higher microdiversity than their virulent counterparts in eight of 12 datasets, whereas the remaining four datasets did not show a significant trend. The detected pattern holds across multiple quality thresholds and lifestyle prediction methods. These findings suggest that lysogeny may promote or preserve genetic variation within phage populations, with potential implications for phage-host coevolution and environmental adaptability.